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Ch ' Abstract 

We consider construction of empirical formulas for predicting new av- 
^ , eraged experimental data from a finite number of the old ones so as to 

■ point out the relevant technical problems. The main problem is that there 

is always only a finite number of inaccurate data at our disposal. As an 
example, we construct empirical formulas for describing Planck's law in 
^ C"j terms of averaged experimental data. 
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1 Introduction 



Physics attempts to predict results of new measurements on the basis of already 
obtained experimental results. So to say, it is concerned with interpolation and 
extrapolation of experimental data by mathematical formulas, starting with 
empirical ones. Historically, many empirical relationships have been stepping 
stones to theories providing physical laws that generalize and extend them. Let 
us give just four examples: 

1) Ptolemaic system espoused in the 2nd century AD, evolved over Coperni- 
can system to Newton's law of universal gravitation, cf. e.g.p]. 

2) Hooke's law stated in 1676 has been generalized to a tensor expression to 
study deformations of various materials, cf. e.g. [2J. 

3) Ohm's law published in 1827 remains an extremely useful formula in elec- 
trical/electronic engineering, cf. e.g. [5] 

4) Balmer formula proposed in 1855 to predict the spectral line emissions of 
the hydrogen atom, and its generalization, Rydberg's formula proposed in 
1888 to predict spectral lines of hydrogen-like atoms of chemical elements, 
were incorporated into the Bohr's model in 1913; a stepping stone to 
quantum mechanics, cf. e.g. [HE] 

How empirical observations of physical behaviour and empirical laws determine 
the form and content of physics and its theoretical structurs has been considered 
in straightforward untechnical style by Cook [5] . Dealing with basic aspects of 
the general framework of physics, he provides a lucid examination of issues 
of fundamental importance in the empirical and metaphysical foundations of 
physics. As a supplement, we will enumerate only the activities necessary when 
constructing empirical formulas based on experimental data. We briefly describe 
the kind and amount of work that is necessary to this end. There are two crucial 
facts: there is always only a finite number of experimental data, and they are 
not accurate. To avoid uniluminating complications, we will consider only the 
averaged experimental data. 

To point out the problems we encounter using experimental data to infer 
empirical laws, let us simulate an empirical derivation of Planck's law. This law 
describes the spectral radiance of electromagnetic radiation of a black body as 
a function of frequency v and its temperature T: 



Take a set of J distinct blackbody states, with Tj denoting the temperature 
of state j = 1, . . . , J. Let us observe the spectral radiance of each state at 
I different frequencies Ui, i = 1, . . . , 7, repeating each measurement N times, 
cf.[6], Sec. 1.2. We denote the actually measured spectral radiance by ij(vi;ri), 




(1) 
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n = 1, . . . , N, and their averages by 



N 

i j (v i )=N- 1 J2 i j(^n)- ( 2 ) 

n=l 

Having obtained a finite number of averaged experimental data about the 
spectral radiance of the j-th blackbody state, say 

D j = {i J (y i ),i = l t ...,l} ) (3) 

we are going to consider two problems: 

(A) How to construct empirical formulas by using only a few of the averaged 
experimental data ij{v%) to predict the rest of them. 

(B) How to construct empirical equations that represent the physical law that 
underlies Dj] that is, how to obtain some quantitative information about 
Planck's law from experimental data. 



2 Empirical formulas 

An empirical formula is a mathematical equation whose parameters have been 
calculated from a few experimental data so as to predict the rest of them. Such 
empirical formulas should be constructed parsimoniously [8], trying to achieve 
the required accuracj0 of predictions using as few parameters as possible. 
Experience suggests we start with a linear ansatz 

Ei(v) = Cl +c 2 y. (4) 

There are many ways to specify parameters c± and ci in terms of ij{yi). Fitting 
two of them, say ijiyi) and ^ (^2), we get 

Ei(v) = L{v\ij(v),vi,vz) , (5) 

where 

L{y;ij(v), v\, i> 2 ) = + [^(^2) - - v\)j{vi - v\) ■ (6) 

This empirical formula reproduces ij(v) well if v is sufficiently close to v\ or 
When a linear ansatz does not reproduce the averaged experimental data 
well enough, we might try a higher-order polynomial ansatz with more parame- 
ters. However, if too many parameters are required, different empirical formulas 
might be more appropriate. 

1 This accuracy should be chosen with regard to the accuracy of measured spectral radiances 
ij(ui;n). 

2 For simplicity, and to avoid being tedious, we will not formally specify what is meant by 
these words. 
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By plotting Iiiia ij{vi) versus Vi, we could conclude that 

Eoo(v) = v 3 exp[L(is; lm/~ 3 ij(p), v u v 2 )\ (7) 

is a suitable^ empirical formula for predicting the asymptotic values of ij {v) for 
large frequencies v >— * oo if we choose v\ and v 2 large enough. 2 

Similarly, on plotting hxijiyi) versus lnz^, we could conclude that the em- 
pirical formula 

Eq(v) =i 3 (is )(is/is ) 2 (8) 

is suitable for predicting the values of ij{vi) for small frequencies Vi — » if we 
choose vq small enough. 2 

One may combine empirical formulas to obtain a more general and parsimo- 
nious empirical formula for predicting the averaged experimental data. 



3 Local empirical approximations of Planck's law 
for radiance of a blackbody state 

As the preceding three empirical formulas turn out useful for predicting the 
averaged experimental data Dj about the radiance of the j-th state, we con- 
struct related local approximations to the underlying physical law. No such 
local approximation can be proved or falsified by the results of a finite number 
of available experimental measurement. They can only show how compatible it 
is with them. 

We start with two basic assumptions about the measurement uncertainties 
of data Dj-. 

Al The averaged experimental data ij(wj) tend toward the values of spectral 
radiance as the number of repeated measurements N becomes very large, 
i.e., for any i, j 

lim ij(ui) = IivuTj) e (0,oo). (9) 

N— »oo 

A2 The maximal relative observational error of data Dj tends to zero as the 
number of repeated observations N becomes large, i.e., 

max|i i (^)-% i ,T ; ,-)|/%i,T,)-*-Q as N -> oo . (10) 

id 

In addition, we make the following three assumptions about the limiting be- 
haviour of the spectral radiance I{y,T): 

A3 The limit 

cij = lim [I(u 2 ,Tj) - I(ui,Tj)]/(v 2 - fi) (11) 

V2—*Vl 

exists for all v\ > 0. 

3 There is a computationally efficient method |51 ll(J| that generalizes such graphical proce- 
dures. Using it, we can also discover more intricate functional dependencies than those that 
can be discovered using graph paper. 
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A4 There is the limit 

C2i = VnI(v,T j )y. (12) 

A5 There is the limit 

Cooj = - lim v~ x ln[/(i/,T,)/i/ 3 ] . (13) 

A6 There is the limit 

c 3j = lim u^e c<xliU ij{u) . (14) 

V — >OG 

These limits define four constants that describe tha intensive physical proper- 
ties of the j-th blackbody state. Evaluation of these limits requires an infinite 
number of averaged experimental data, but the number of available experimen- 
tal data is always finite. Therefore, such limits are theoretical constructs that 
are evaluated by how useful they are in making empirical approximations to 
the underlying physical law. If, in some application, we need the actual value 
of such a limit, we have to make do with an estimate based on the available 
experimental data Dj, and denote it with the same letter written with a hat: 
cij, &2j, C3j, dooj. Metrological methods provide various estimates to this end. 

Inspecting the values of c 3 j, and of the product C2jC OCJ j for the states consid- 
ered, we might conclude that the values c 3 j and C2jC oa j are state- independent 
and equal the same value, say c 3 , i.e. that 

C2jC ooj = c 3j = c 3 for all j = 1, 2, . . . , J . (15) 

3.1 Linear approximation 

The empirical formula E\{u) suggests the following linear approximation to 
Planck's law: 

Li(v,Tj) = I{u u Tj) + cij-(i/-i/i) + 0{{y-v 1 f) as v -f v x . (16) 

Testing it on the available data Dj, we get an estimate c\j of the value of 
constant cij, and an estimate of for how large values of \v — v\\ this linear 
approximation is still useful. Hooke's law and Ohm's law are two examples of 
such linear approximations. 

3.2 Asymptotic behaviour for low frequencies 

The uscfullncss of the empirical formula E n {v) indicates that spectral radiance 
I(u, Tj) has the following asymptotic behaviour for low frequencies: 

L (v,Tj) = c 2 jv 2 + 0(u 3 ) as ^^0. (17) 

The first term, the Rayleigh-Jeans formula, is totally useless for predicting the 
spectral radiance for high frequencies. 
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3.3 Asymptotic behaviour for high frequencies 

The uscfullncss of the empirical formula E^iy) for predicting the values of ij(vi) 
for high frequencies suggests that 

L O0 (v,T 3 ) = c 3 v 3 e- c ^ 1/ + 0(v 3 e- 2c °°^) as v -► oo (18) 

describes the asymptotic behaviour of I(i/,Tj) for high frequencies. The first 
term, Wien's law, predicts also the correct limiting value of spectral radiance for 
v = 0, but does not imply the correct asymptotic behaviour for low frequencies. 

4 Global approximations to Planck's law for the 
radiance of a blackbody state 

As the local approximations Li, L and turn out to predict the averaged 
experimental data Dj well enough locally, we will use them to construct a global 
approximation to Planck's law for the j-th state. They suggest we assume that 

A7 Spectral radiance is a non-negative, analytic function in a vicinity of any 
positive frequency. 

After some juggling, we construct a mathematical formula that unifies local, 
asymptotic approximations L and L m by an interpolation and agrees with 
assumption A7: 

Ufa C3, Cooj, Cj) = c 3 v 2 {cj + v) j [e c ^ v - 1 + C jCooj ) , (19) 

where Cj is a non-negative empirical constant. On testing this formula with 
Cj = 0, we find that it represents data Dj well enough for any state if N is 
large enough. So we could presume that U (y \ c 3 , Cooj , 0) represents Planck's law 
for the radiance of the j-th blackbody state; it is determined by two physical 
properties c 3 and 0^ of this state. 

5 Empirical form of Planck's law 

On measuring the temperatures Tj of various states and plotting TJ 1 versus 
Cooj, we might conclude that there is a positive constant Coo such that 

c 00j =c 00 T7 1 ,j = l,2,...,J. (20) 

And so we could put forward the formula 

U(v;c 3 , Cqc /T,0) (21) 

as an empirical form of Planck's law that is valid for all temperatures T and 
depends only on two empirical constants, c 3 and c^. 
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6 Comments 



6.1 Identifying the states 

The blackbody states are identified by the index j, giving the order in which 
they were researched. In principle, we could identify the blackbody states by 
their physical properties. So, were they known to us, we could use the values 
of Cooj that are inversly proportional to temperature, or the values of spectral 
radiance for some frequency v m . However, we could make do also with estimates 
Cooj and I{v m , Tj), were they so accurate that they did not overlap for the states 
considered. 

6.2 Qualitative properties of Planck's law 

Planck's law is such that 

I(-v,-T) = -I{v,T), I(Tv,T)=T 3 I(v,l), (22) 

and I(y, T) is an analytic function of the complex variable v everywhere but at 
hv = 2mrkT 7 n = 1, 2, . . . , where it has a first-order pole, and at \ v\ — oo, where 
it has an essential singularity. However, these properties cannot be directly 
inferred from experimental data. 

Note that only the constant function is analytic in the whole complex plane. 
So, an analytic funtion is, up to a constant, uniquely determined by its singular- 
ities. An estimate of a singular point of l(y, T) close to a point on the positive 
real axis can be obtained from the convergence of the Taylor expansion at this 
point. 

Were some physical law an even function of v, this property would not be di- 
rectly evident from experimental data, though equations that presume this prop- 
erty would be more efficient, especially so in the vicinity of v = 0. The asymp- 
totic behaviour L (v,Tj) is compatible with the hypothesis that Planck's law 
is an even function of frequency. It is only the asymptotic behaviour L^v, Tj) 
that falsifies this hypothesis. 

6.3 The time we have for experiments is limited 

The maximal relative observational error of the available data Dj is in general 
positive, though we expect it to get smaller as we increase the number N of 
repetitions. To obtain J sets Dj, each containing / data, we must perform NIJ 
mesurements. And each one takes a certain amount of time, say at least t m . 
However, there is obviously an upper limit, say Tm, on the amount of time 
available for observing physical phenomena. So the total number of possible 
measurements is at most TM/t m , and we must have 

UN < T M /t m . (23) 

Let us point out some consequences of this limit on the possible number of 
available experimental data. 
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A) When T m < t m , the phenomenon is practically unobservable; an example 
is provided by the blackbody radiation of visual frequencies at room tempera- 
tures. So parts of Planck's law are practically untestable, cf. e.g. [7]. 

B) For each physical law L(x), there is an infinite number of alternatives 
L a (x) = L(x)(l + e(j)(x)), with 4>(x) bounded: so the relative difference 

\(L(x) - L a (x))/L(x)\ < |e| sup \(f>(x)\ (24) 

X 

is arbitrarily small for sufficiently small |e|. And when relative differences be- 
tween predictions of two physical laws are sufficiently small, we will never be 
able to tell them apart experimentally. Such a problem occurs with the empirical 
form of Planck's law (|2Tj) , since the relative difference between U{v : C3, Coo/T, Cj) 
and /7(^; C3, Coo/T, 0) tends to zero uniformly as Cj — ► 0. Thus we will never 
be able to falsify the assumption that Cj is a very small positive constant. It is 
expedient to choose Cj = 00 though U{v; C3, Cqo/T, Cj ^ 0) might turn out to 
be a theoreticaly significant modification of Planck's law. 
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4 Such is the case with the mass of the photon. For some theoretical calculations it is 
convenient to limit the photon mass to zero only in their final stage, though many theoretical 
considerations take as their basic presumption that photons have no mass, cf. Veltman|lf I. 
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